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Abstract 

While the X-ray, GeV gamma-ray, and TeV gamma-ray skies have been extensively studied, the MeV gamma-ray 
sky is not well investigated after the Imaging Compton Telescope (COMPTEL) scanned the sky about two decades 
ago. In this paper, we investigate prospects for active galactic nuclei population studies with future MeV gamma- 
ray missions using recent spectral models and luminosity functions of Seyferts and flat spectrum radio quasars 
(FSRQs). Both of them are plausible candidates as the origins of the cosmic MeV gamma-ray background. If the 
cosmic MeV gamma-ray background radiation is dominated by non-thermal emission from Seyferts, the sensitivity 
of 10 -12 erg cm -2 s _1 is required to detect several hundred Seyferts in the entire sky. If FSRQs make up the 
cosmic MeV gamma-ray background, the sensitivity of ~ 4 x 10“ 12 erg cm -2 s -1 is required to detect several 
hundred FSRQs following the recent FSRQ X-ray luminosity function. However, based on the latest FSRQ gamma- 
ray luminosity function, with which FSRQs can explain up to ~ 30% of the MeV background, we can expect several 
hundred FSRQs even with the sensitivity of 10” 11 erg cm -2 s _1 which is almost the same as the sensitivity goal of 
the next generation MeV telescopes. 
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1. Introduction 

The X-ray, GeV gamma-ray, and TeV gamma-ray skies have 
been deeply investigated by various telescopes. However, the 
MeV gamma-ray sky has not been extensively studied yet af¬ 
ter the Imaging Compton Telescope (COMPTEL) onboard the 
Compton Gamma-Ray Observatory ( CGRO ) satellite scanned 
the sky in 1990s (Schonfelder et ah, 2000). This means that 
huge discovery space is buried in this energy range. 

Measurement of MeV gamma rays is hampered by huge 
background events. Under this difficulty, the COMPTEL 
achieved the sensitivity of ~ 10~ 10 erg cm~ 2 s _1 in the MeV 
band and discovered ~63 gamma-ray sources (32 of them 
were gamma-ray bursts) at 0.75-30 MeV (Schonfelder et ah, 
2000). The COMPTEL has found 10 extragalactic active 
galactic nuclei (AGNs) and 5 unidentified high latitude objects 
(Schonfelder et ah, 2000). 

To open a new observational window in the MeV sky, vari¬ 
ous missions have been proposed. ASTRO-H (Takahashi et ah, 
2014), which will be launched in 2015, will have a sub-MeV 
instrument, the soft gamma-ray detector (SGD, Tajima et ah, 
2010; Watanabe et ah, 2012) covering at 40-600 keV. The 
SGD will provide precise spectra with higher sensitivity than 
previous instruments by introducing the concept of narrow- 
Field-of-View Compton telescope which can dramatically re¬ 
duce background events (Takahashi et ah, 2001, 2004). Other 
instruments are also proposed such as CAST (Nakazawa et ah, 


2012), DUAL (von Ballmoos et ah, 2012), GRIPS (Greiner 
et ah, 2012), SMILE (Takada et ah, 2011), COSI (Kierans 
et ah, 2014), and ASTROGAM 1 . Moreover, the instrumen¬ 
tal performance were tested with ballon experiments (Takada 
et ah, 2011; Bandstra et ah, 2011). These future missions will 
unveil the MeV sky in near future. In addition, a new image 
reconstruction algorithm that has improved sensitivity to mul¬ 
tiple point-like sources has been proposed toward those future 
all-sky missions (Ikeda et ah, 2014). The new era of the MeV 
gamma-ray astronomy is approaching. 

A fundamental question to those missions is how many 
sources one can expect with their designed instruments. In the 
extragalactic sky, the dominant population in the MeV sky is 
still uncertain because the origin of the cosmic MeV gamma- 
ray background from 0.3 MeV to 100 MeV is not well estab¬ 
lished (see Inoue, 2014, for a recent review). Seyferts which 
make up the cosmic X-ray background (CXB) cannot explain 
this MeV background because of the assumed cutoff at a few 
hundred keV (e.g. Ueda et ah, 2003; Gilli et ah, 2007). Above 
100 MeV, the total gamma-ray background is known to be pre¬ 
dominantly composed of blazars (e.g. Inoue & Totani, 2009; 
Ajello et ah, 2015), star-forming galaxies (e.g. Ackermann 
et ah, 2012), and radio galaxies (e.g. Inoue, 2011). Although 
the MeV background may comprises these GeV gamma-ray 
populations as well, the lower energy part of the MeV back- 

1 http://astrogam.iaps.inaf.it/index.html 
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ground spectrum is smoothly connected to the CXB spectrum 
and shows softer (photon index of T ~ 2.8 Fukada et al., 1975; 
Watanabe et ah, 1997; Weidenspointner et ah, 2000) than the 
GeV component (photon index of F ~ 2.4 Ackermann et ah, 
2015a), indicating their different origins. 

Several candidates have been considered as the origin of the 
MeV background. One was the nuclear-decay gamma rays 
from Type la supernovae (Clayton & Ward, 1975; Zdziarski, 
1996; Watanabe et ah, 1999). However, the measurements 
of the cosmic Type la supernovae rates show that the rates 
are not enough to explain the observed flux (e.g. Ahn et ah, 
2005; Strigari et ah, 2005; Horiuchi & Beacom, 2010; Ruiz- 
Lapuente et ah, 2015). Seyferts may naturally explain the 
MeV background up to a few tens of MeV and the smooth 
connection to the CXB (Schoenfelder, 1978; Field & Rogers, 
1993; Stecker et ah, 1999; Inoue et ah, 2008, hereinafter 108). 
Non-thermal electrons in coronae can generate the power-law 
spectrum in the MeV band after the thermal cut-off via the 
Compton scattering of disk photons (108). Flat-spectrum ra¬ 
dio quasars (FSRQs) whose peak in the spectrum locates at 
MeV energies (Blom et ah, 1995; Sambruna et ah, 2006) are 
also expected to significantly dominate the MeV background 
(Ajello et ah, 2009, hereinafter A09). Radio galaxies have been 
also discussed as one of the origins of the MeV background 
(Strong et ah, 1976). Recent studies have revealed that both of 
lobe (Massaro & Ajello, 2011) and core (Inoue, 2011) emis¬ 
sions from radio galaxies could contribute only ~ 10% of the 
MeV gamma-ray background flux. Dark matters has also been 
discussed as the origin of the MeV background (e.g. Ahn & 
Komatsu, 2005a,b). Those MeV mass dark matter particle can¬ 
didates are less natural than GeV-TeV dark matter candidates. 
Therefore, Seyferts and FSRQs can be regarded as potential 
astrophysical origins of the MeV background. 

The purpose of this paper is to investigate prospects of ex- 
tragalactic observations by future MeV instruments, especially 
of statistical aspects of AGNs. We focus only on Seyferts and 
FSRQs in this paper. For this purpose, we use recent luminos¬ 
ity functions (LFs) and spectral models of Seyferts and FSRQs 
in literature. Since a quantitative estimate based only on the 
COMPTEL results is not easy due to the small number of de¬ 
tected sources, we adopt recent X-ray or gamma-ray luminosity 
functions and spectral models of Seyferts and FSRQs. 

This paper is organized as follows. MeV gamma-ray emis¬ 
sion models of Seyferts and FSRQs are described In Section 
2. LFs of Seyferts and FSRQs are presented in Section 3. In 
Section 4, results of the expected source counts and redshift 
distributions are presented. Discussions and conclusions are 
given in Section 5. Throughout this paper, we adopt the stan¬ 
dard cosmological parameters of (/i,Qm,^a) = (0.7,0.3,0.7). 

2. Spectra of Active Galactic Nuclei in the MeV gamma- 
ray band 

2.1. Seyferts 

The X-ray spectra of Seyfert are phenomenologically ex¬ 
plained by a combination of following components: a primary 
power-law continuum with a cutoff at ~0.3 MeV in the form 
of J5 _r exp(— E/E c ), absorption by surrounding gas, emission 
lines, a reflection component, and a soft excess of emission at 
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Fig. 1 . Spectral templates for Seyferts and FSRQs. Solid, dashed, dot¬ 
ted, and dot-dashed curve corresponds Ueda et al. (2003, U03) Seyfert 
spectral model with a thermal cutoff at 0.3 MeV with Njj = 10 20 cm 2 , 

U03 Seyfert spectral model with a thermal cutoff at 0.3 MeV with 
Ah = 10 25 cm 2 , Inoue et al. (2008,108) Seyfert spectral model with 
a non-thermal tail in the MeV band with Njj = 10 20 cm 2 , and Ajello 
et al. (2009, A09) FSRQ spectral model having a peak at 1 MeV. The 
shape is assumed to be independent of luminosities for these models. 

)$ 2 keV (e.g. Dadina, 2008). Relative fractions of these com¬ 
ponents vary with sources. Physically Comptonization of disk 
photons in a corona above the accretion disk generate the pri¬ 
mary power-law continuum (see e.g. Katz, 1976; Pozdniakov 
et al., 1977; Sunyaev & Titarchuk, 1980). The temperature of 
the corona roughly determines the position of the spectral cut¬ 
off and the photon index of the intrinsic continuum together 
with the optical depth (see e.g. Zdziarski et al., 1994). The 
Compton reprocessed emission and bound-free absorption of 
the primary continuum by surrounding cold matter generate the 
reflection component (Lightman & White, 1988; Magdziarz & 
Zdziarski, 1995). 

Although coronae with temperature of ~ 0.1 MeV do not 
produce significant MeV gamma-ray emission, some spectral 
models predict a power-law tail after the thermal cutoffs (e.g. 
108). If the corona is composed of thermal and non-thermal 
populations, a MeV tail will appear after the cut-off. Such non- 
thermal electrons may exist if magnetic reconnection heats the 
corona (Liu et al., 2002) as in the Solar flares (e.g. Shibata 
et al., 1995) and Earth’s magnetotail (Lin et al., 2005). The 
MeV background can be explained by the Seyferts explaining 
the CXB which have non-thermal electrons in the coronae hav¬ 
ing ~ 4% of the total electron energy (108). 

Observationally, the cutoff energy of Seyferts are con¬ 
strained at > 0.2 MeV (Ricci et al., 2011). The Oriented 
Scintillation Spectroscopy Experiment ( OSSE ) clearly detected 
emission up to 0.5 MeV in the spectrum of the brightest Seyfert 
NGC 4151 (Johnson et al., 1997). This measurement con¬ 
strained non-thermal fraction to be <15% (Johnson et al., 
1997). Interestingly, future radio observations are able to probe 
those non-thermal tail through synchrotron emission (Laor & 
Behar, 2008; Inoue & Doi, 2014). 

In this paper, as shown in Figure 1, two primary spectral 
models are applied for Seyferts as in Inoue et al. (2013). One 
is thermal spectral model having a power-law continuum with 
a cutoff (see e.g. Ueda et al., 2003, hereinafter U03). We adopt 
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Fig. 2. Spectral templates for FSRQs, but for A12 models. 
Solid, dashed, dotted, and dot-dashed curve corresponds to the 
0.1-100 GeV luminosity of L-y = 10 48 ' 9 , 10 48 ' 2 , 10 47 ’ 8 , and 
10 46 ' 8 ergs -1 , respectively, based on Ajello et al. (2012). 

r = 1.9 and E c = 0.3 MeV. The other is thermal plus non- 
thermal spectral model (see 108 for details) which could ex¬ 
plain the MeV background spectrum. We adopt the same pa¬ 
rameters as in 108 with the non-thermal photon index of 2 . 8 , 
but setting the cutoff at 0.3 MeV. For the Compton reflection 
component, we use a Compton reflection model (Magdziarz & 
Zdziarski, 1995), assuming a solid angle of 27 t, an inclination 
angle of cos i = 0.5, and solar abundance for all elements. To 
calculate absorbed spectra, we use an absorption model called 
”wabs” developed for the XSPEC package. 

2.2. FSRQs 

Multi-wavelength observations allowed us to study overall 
spectra of blazars from radio to gamma-ray. The observed 
blazar spectra consists of two broadband emission components 
(Ulrich et al., 1997; Fossati et al., 1998; Ghisellini et al., 1998; 
Kubo et al., 1998; Abdo et al., 2010). One is synchrotron ra¬ 
diation component and the other is inverse Compton (IC) com¬ 
ponent in which electrons scatter internal synchrotron radia¬ 
tion (see e.g. Jones et al., 1974) or external radiation (see e.g. 
Dermer & Schlickeiser, 1993; Sikora et al., 1994). 

Blazars are divided into two categories by their optical spec¬ 
tra. Those are BL Lacertae objects (BL Lacs) and FSRQs. In 
the case of BL Lacs, the highest observable energy part of the 
synchrotron emission appears in the X-ray band. X-ray spectra 
of BL Lacs show soft spectra with a photon index of T ~ 2 — 3. 
In the case of FSRQs, the lowest observable energy part of the 
IC emission appears in the X-ray band showing harder spectra. 
Although BL Lacs are not expected to the dominant contribu¬ 
tors of the MeV background, FSRQs may make up the whole 
MeV background (A09). Hence, in this paper, we do not con¬ 
sider BL Lacs but only FSRQs for blazars. 

We consider two spectral energy distribution (SED) mod¬ 
els for LSRQs. One is based on the Swift-BAA FSRQs (A09), 
while the other is based on the Fermi -LAT FSRQs Ajello et al. 
(2012, hereinafter A12). Following A09, the SED is given by 
a double power-law model: 



(1) 


where we set Eb = 1 MeV, Ti = 1.6, and T 2 = 2.9 following 
A09 (see Figure 1). The average photon index of the Swift- 
BAT FSRQs is 1.6 ±0.3 at 15-55 keV (A09). Although the¬ 
oretical models predict a spectral turnover at the MeV band 
for FSRQs (see e.g. Inoue & Takahara, 1996), Eb and L ’2 in 
the MeV band are observationally unconstrained. Hence, these 
parameters are set FSRQs to make up the MeV background. 
Here, the latest catalog of the Fermi -LAT sources reported de¬ 
tection of 467 FSRQs. The mean of their photon index above 
0.1 GeV of those gamma-ray bright FSRQs is typically ~ 2.4 
(Ackermann et al., 2015b), which is harder than that of the 
MeV background. Combining with the Fermi -LAT and Swift- 
BAT detected FSRQs, the average SED is modelled in Fig. 9 of 
A12 which are shown in Figure 2. We also adopt these SEDs 
as the other FSRQ template. However, as discussed in the next 
section, the whole MeV background is not be able to be ex¬ 
plained by the A12 FSRQ model. 


3. Luminosity Functions of Active Galactic Nuclei 


To obtain the expected source counts and redshift distribu¬ 
tions of AGNs in the MeV gamma-ray band, a LF is required. 
We briefly review LFs of Seyferts and FSRQs in this section. 

3.1. Seyferts 

The cosmological evolution of Seyferts is investigated by 
various X-ray surveys (see e.g. U03; Hasinger et al., 2005; 
Gilli et al., 2007; Ueda et al., 2014). Based on these surveys, 
the observed X-ray luminosity functions (XLFs) are known to 
be well represented by luminosity-dependent density evolution 
(LDDE) models (see e.g. Ueda et al., 2014, and references 
therein). 

In this study, we follow the U03 LDDE XLL at 2-10 keV 
for the Seyfert models, since the 108 Seyfert model is based on 
U03. We note that U03 modelled the evolution of X-ray emit¬ 
ting AGNs. The comoving number density p x in the LDDE 
is 


Px(Lx,z,Nk) = p x (L x ,0) f (L x , z)ifNbi', L x , z), (2) 


where L x is the X-ray luminosity at 2-10 keV before the ab¬ 
sorption, z is the redshift, and IVh is the neutral hydrogen col¬ 
umn density. p x (L x , 0) is the AGN XLL at present day char¬ 
acterized by the faint-end slope index 71 , the bright-end slope 
index 72 , and the break luminosity L*, as: 


Px(L x ,0) = A 




(3) 


where A is the normalization parameter having a dimension of 
volume -1 . 

The density evolution function f(L x ,z) is given as: 


f(L x ,z) = 


(1 + z) Pl z<z c (L x ), 

(l + Zc (L x )r{ TT ^) P2 z>z c {Lx\ A) 


where z c is the redshift of evolutionary peak, given as 
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Table 1. The parameters of the AGN Luminosity Functions 



Ueda et al. (2003) 
Seyfert 

2-10 keV 

Ajello et al. (2009) 
FSRQ 

15-55 keV 

Ajello et al. (2012) 
FSRQ 

0.1-100 GeV 

A a 

50.4 ±3.3 

0.533 ±0.104 

0.0306 ±0.0023 

l°gic)7 j * fc 

43.94l°; 24 

44.0 

log 10 (0.84±0.49)+48.0 

7i 

0.86 ±0.15 

- 

0.21 ± 0.12 

72 

2.23 ±0.13 

3.45 ±0.20 

1.58 ±0.27 

* 

z c 

1.9 

- 

1.47 ±0.16 

logio La 

44.6 

- 

- 

a 

0.335 ±0.07 

- 

0.21 ±0.03 

pi 

4.23 ±0.39 

3.72 ±0.50 

7.35 ±1.74 

P2 

-1.5 

-0.32 ±0.08 

-6.51 ±1.97 

P 

- 

- 

2.44 ±0.01 

a 

- 

- 

0.18 ± 0.01 


a : In units of 10 7 Mpc 3 . 
b : In units of erg/s. 


z c (Lx) = 


z*(L x /L a y 


Lx > L a , 
Lx < L a . 


(5) 


The function ^(TVh; Lx, z) is the distribution of absorption 
column density given in the following form in the XLF (Ueda 
et al., 2003): 

2 -^y(Lx,z) 20.0 < log Ah < 20 


77(IVh; Lx,z)={ 1 +zi>(Lx, z) 

I ^ip(Lx,z) 


20.5 < logA H < 23 
23.0 < log Ah < 24 


where e =1.7 and t(>(Lx, z) = min{t/j max , max[0.47 — 
0.1(logL x -44.0),0])}, for which t/Wx= (1 + e)/(3 + e). We 
note the integration of t/(Ah;L x ,z) between 20.0 < log Ah < 
24.0 is unity. 

The parameters obtained by the fit to the observed data of X- 
ray AGNs in U03 are shown in Table 1. We set the minimum 
and maximum of the X-ray luminosity as Lx, min = 10 415 and 
Lx, max = 10 48 ° erg s^ 1 , respectively, the same as in U03. 
The redshift range is set to be 0.0 < z < 5.0, otherwise noted. 
To explain the CXB, the fraction of the Compton thick AGNs 
between 24.0 < log Ah < 25.0 is set to be the same as that of 
the population at 23.0 < log An < 24.0 following U03. 

Recently, Aird et al. (2010) suggested another evolution 
model for AGNs, the luminosity and density evolution (LADE) 
model. Although the distribution of absorption column den¬ 
sity is not available for the LADE, we can test that model 
assuming the same absorption column density distribution as 
in U03. The overall source counts in the MeV band signifi¬ 
cantly decrease by a factor of ~ 3 at the flux range of 10” 11 - 
10 -12 erg cm -2 s _1 . However, the latest study of X-ray AGN 
evolution (Ueda et al., 2014) found that the AGN XLF is not 
well described with the LADE model but with LDDE model. 
When we adopt the latest LDDE model (Ueda et al., 2014), the 
Seyfert source counts in the MeV band decrease by a factor of 
~ 40% at the flux range of 10 _11 -10 -12 erg cm -2 s _1 . As the 
108 Seyfert model is based on U03, we follow the U03 LDDE 
model in this paper. 


3.2. FSRQs 

In the GeV gamma-ray band, A12 and Ajello et al. (2014) 
have recently shown that the gamma-ray luminosity functions 
(GLFs) of FSRQs and BL Lacs can be described by LDDE 
models based on the Fermi -LAT FSRQs and BL Lacs, which 
gwas indicated before the launch of Fermi (e.g. Narumoto & 
Motani, 2006; Inoue & Totani, 2009). With the A12 GLF 
ythodel, FSRQs can explain up to ~ 30% of the MeV back¬ 
ground. In the X-ray band, A09 have investigated the XLF of 
FSRQs using the 3-year Swift-BAT AGN survey data including 
26 FSRQs. A pure luminosity evolution (PLE) model success¬ 
fully reproduced the observed distribution of X-ray FSRQs. 
Assuming this PLE model and the spectral model in §. 2.2, 
FSRQs explained the whole MeV background. 

These two recent FSRQ XLF and GLF models of A09 and 
A12 are adopted in this paper. The FSRQ XLF at 15-55 keV 
by A09 is given in the form of the PLE. The comoving number 
densities px in the PLE is given as: 


Px{L x ,z) = px(L x /e(z), 0), 


(7) 


where Lx is the X-ray luminosity at 15-55 keV here. The local 
XLF is described as: 


Px(L x ,0) = A 


Lx 

T* 


( 8 ) 


The evolution factor e(z) is independent of luminosity given 
by: 

e(z) = (1 + z) Pl+P2Z (9) 

The values of the parameters for the A09 XLF are shown in 
Table 1. The minimum luminosity is also assumed to evolve 
with redshift as 


Lx,min(z) — Lx, to. i n ,0 X e(z), 


( 10 ) 


where Lx,min, o is the minimum luminosity at z = 0. We set 
Lx,min,o = 3 x 10 44 erg s _1 , Lx, ma x=10 50 0 erg s _1 , 0 min = 
0.0, and z max = 6.0 (A09), otherwise noted. 

On the contrary, the A12 FSRQ GLF is given in the form 
of LDDE model. The comoving number densities p 1 in the 
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U03 Seyfert 

108 Seyfert - 

A09FSRQ 
A12FSRQ - 


Fig. 3. Cumulative source counts in the entire sky as a function of MeV gamma-ray flux (in S ) of Seyferts and FSRQs. The five pan¬ 
els correspond to different photon energy bins as indicated in the panels. The solid, dashed, dotted, and dot-dashed curves are predic¬ 
tions based on Ueda et al. (2003, U03) for Seyferts, Inoue et al. (2008, 108) for Seyferts, Ajello et al. (2009, A09) for FSRQs, and Ajello 
et al. (2012, A12) for FSRQs. The horizontal thin double-dot-dashed line corresponds to one source detection level. Due to the ther¬ 
mal cutoff at 0.3 MeV in the spectral model of U03, detection of Seyferts above 3 MeV is not expected in the given flux ranges. 


LDDE is given as: 

p- y (L 1 ,z,T)=p 1 (L^,0)f(L y ,z)e(T), (11) 


where L 7 is the gamma-ray luminosity at 0.1-100 GeV and I’ 
is the intrinsic photon index between 0.1-100 GeV. p 7 (L 7 ,0) 
is the FSRQ GLF at present. This is characterized as: 


p 7 (i 7 ,0) — A 




( 12 ) 


The density evolution function /(i 7 , z) is given by: 


f(Ly,z) 


1 + z 

1 + z c (L 7 ) 



1 + z 

1 + z c (Lj) 



where 2 C is given as z c {L 1 ) = z* ■ (L 7 /10 48 )“ 

The function 0(F) describes the distribution of photon index, 
which is given by a Gaussian form in the GLF (A12): 


0 ( 7 ) = exp 


(T -M) 2 ' 

2 a 2 


(14) 


where fi and cr are the mean and the dispersion of the Gaussian 
distribution. However, we adopt the average SED of FSRQs 
(see Fig. 2) not power-law SEDs. Therefore, we include the 
integrated 0(T) over T between 1.8 and 3.0 which is indepen¬ 
dent of luminosity and redshift. 

The parameters for the A12 GLF are shown in Table 1. The 
limits of gamma-ray luminosities, redshifts, and photon in¬ 
dices are set as L 7)ln ; n =10 44 erg s _1 , L 7 , m i n =10 52 erg s _1 , 
Zmin=0.01, and z m ax=6, the same as in A12, otherwise noted. 


4. Source Counts and Redshift Distribution 

The cumulative source counts in the entire sky at a sensitiv¬ 
ity of S, N(> S), are given by integrating the luminosity func¬ 
tion where the luminosity range is set to be L ul [ n (SE 1 ,E 2 +) < 
L < L max . Se 1 ,e 2 represents the sensitivity in a given energy 
range between Ei and 1+ in the unit of [erg cm -2 s' 1 ] and 
L m in{SE 1) E 2 , z) is the minimum detectable luminosity with 
a given sensitivity SE lt E 2 and a redshift 2 . Figure. 3 shows 
the cumulative source counts in the entire sky as a function 
of sensitivity limit at a given energy band in the entire sky 
predicted based on U03 for Seyferts with the thermal cutoff, 
108 for Seyferts with the non-thermal tail, A09 for Swift-BXY 
FSRQs, and A12 for Fenni-LAT FSRQs. The surface number 
density is given by dividing the number by 47 t steradian. The 
energy bands are divided into five as 0.3-1 MeV, 1-3 MeV, 
3-10 MeV, 10-30 MeV, and 30-100 MeV as indicated in pan¬ 
els. Hereinafter, we investigate required sensitivities for detec¬ 
tion of several hundred sources. Several hundred samples may 
allow us to study cosmological evolution of a population (see 
U03 for X-ray AGNs, A12 for gamma-ray FSRQs, Ajello et al., 
2014, for gamma-ray BL Lacs). We also consider the sensitiv¬ 
ity limit of 10 -11 erg cm -2 s^ 1 at the MeV band which is an 
order of magnitude better than the sensitivity of the COMPTEL 
and almost the same as the expected sensitivity of the next gen¬ 
eration MeV instruments such as ASTROGAM. 

Based on the U03 Seyfert model, ~ 20 Seyferts at 0.3- 
1 MeV with the sensitivity limit of 10~ n erg cm -2 s _1 are 
expected, which would be useful to investigate the position of 
the thermal cutoff energies of Seyferts. Once we can achieve 
the sensitivity limit of 10 -12 erg cm' 2 s _1 , ~ 1400 Seyferts 
at 0.3-1 MeV. However, only ~ 10 Seyferts are expected at 1- 
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Fig. 4. Expected cumulative redshift distribution of AGNs in the entire sky assuming the limiting sensitivity of 10“ 11 [erg cm" 2 s“ 1 ]. 
The five panels correspond to different photon energy bins as indicated in the panels. The solid, dashed, dotted, and dot-dashed curves 
are predictions based on Ueda et al. (2003, U03) for Seyferts, Inoue et al. (2008, 108) for Seyferts, Ajello et al. (2009, A09) for 
FSRQs, and Ajello et al. (2012, A12) for FSRQs. The horizontal thin double-dot-dashed line corresponds to one source detection level. 


3 MeV. As the U03 Seyfert model relies on the thermal cutoff 
spectral model at 0.3 MeV, we can not expect any detection of 
U03-type Seyferts above several MeV. 

If the non-thermal tails of Seyferts make up the MeV back¬ 
ground (108), more Seyferts at the MeV band would be ex¬ 
pected. Based on the 108 Seyfert model, we can expect 
~ 60, ~ 10, ~ 5, and a few Seyferts at 0.3-1, 1-3, 3-10, 
and 10-30 MeV, respectively, with the sensitivity limit of 
10“ n erg cm -2 s -1 . By achieving the sensitivity limit of 
10“ 12 erg cm -2 s -1 , we can expect a reasonable number of 
detections as ~ 1800, ~ 400, ~ 140, ~ 40, and ~ 10 Seyferts 
at 0.3-1, 1-3, 3-10, 10-30, and 30-100 MeV, respectively. 
Therefore, the sensitivity limit of ~ 10 “ 12 erg cm -2 s -1 is re¬ 
quired to detect several hundred Seyferts even if Seyferts dom¬ 
inate the MeV background radiation. 

If the origin of the MeV background is FSRQs (A09), we 
can expect ~ 90, ~ 50, ~ 30, ~ 10, and ~ 5 FSRQs at 0.3- 
1, 1-3, 3-10, 10-30, and 30-100 MeV, respectively, with the 
sensitivity limit of 10“ n erg cm -2 s —1 . Once the sensitivity 
limit of 4 x 10 “ 12 erg cm -2 s -1 is achieved, ~ 720, ~ 420, 
~ 260, ~ 90, and ~ 40 FSRQs at 0.3-1,1-3, 3-10, 10-30, and 
30-100 MeV will be detected, respectively, because of strong 
positive evolution. However, based on the latest A12 FSRQ 
GLF, with which the whole MeV background can not be ex¬ 
plained solely by FSRQs, the expected number will be ~ 220, 
- 260, - 390, ~ 370, and ~ 380 FSRQs at 0.3-1, 1-3, 3- 
10, 10-30, and 30-100 MeV, respectively, with the sensitivity 
limit of 10“ n erg cm -2 s -1 . Moreover, at the sensitivity limit 
of the COMPTEL, ~ 10 “ 10 erg cm ' 2 s“\ about 10 FSRQs 
are expected based on the A12 GLF, which is consistent with 
the number of the COMPTEL blazars. The difference between 
FSRQ models mainly comes from the different evolutionary 


history. 

When the sensitivity limit of 10 “ 12 erg cm -2 s -1 is 
achieved in the lowest energy band, expected detection rate for 
Seyferts and FSRQs will be one Seyfert per 22-28 deg 2 and 
one FSRQ per 4-40 deg 2 , respectively, depending on assumed 
models. To avoid the source confusion, we may need angu¬ 
lar resolutions of several degrees. ASTROGAM will achieve 
the FWHM of the angular resolution measure of ~ 4 deg at 
0.3 MeV 2 . We note that COMPTEL had those of ~ 8 deg at 
0.5 MeV (Schoenfelder et al., 1993). 

It is also important to know how distant AGNs can be ob¬ 
served by future missions. Since FSRQs have a spectral peak 
at the MeV band, future MeV missions will be able to detect 
distant FSRQs. In the GeV band, the most distant blazar con¬ 
firmed by Fermi is PKS 0537-286 at z = 3.1 at > 100 MeV 
(Ackermann et al., 2015b). Romani et al. (2004) and Romani 
(2006) have reported one EGRET blazar candidate at z ~ 5.48, 
but it has not been confirmed so far by Fermi. It is also ex¬ 
pected that Fermi may eventually detect blazars at z > 6 (Inoue 
et al., 2011). Takahashi et al. (2013) have reported a candidate 
Fermi gamma-ray blazar at z ~ 3 — 4. 

Figure. 4 shows the cumulative redshift distribution N(> z ) 
as a function of redshift at a given energy band in the entire 
sky assuming the sensitivity limit of 10“ n erg cm ' 2 s _1 pre¬ 
dicted based on U03 Seyfert model, 108 Seyfert model, A09 
FSRQ model, and A12 FSRQ model. Here, we set z ni;ix = 8 as 
the most distant known quasar is at z = 7.085 (Mortlock et al., 
2011). The energy bands were divided into five as 0.3-1 MeV, 
1-3 MeV, 3-10 MeV, 10-30 MeV, and 30-100 MeV as indi¬ 
cated in panels. Figure. 5 shows the same but assuming the 
sensitivity limit of 10“ 12 erg cm -2 s _1 . 

2 http://astrogam.iaps.inaf.it/index.html 
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Fig. 5. The same as Figure. 4, but for the limiting sensitivity of 10 12 [erg cm 2 s 1 ], 


For the Seyfert models, the detectable Seyferts will 
be only local Seyferts. Even if the sensitivity limit of 
10 -12 erg cm -2 s -1 is achieved, the most distant Seyfert 
will be at z ~ 1. On the other hand, for the FSRQ mod¬ 
els, we can expect 2 ~ 3-4 FSRQs with the sensitivity 
limit of 10~ n erg cm -2 s _1 . Once the sensitivity limit of 
10 -12 erg cnr 2 s" 1 is achieved, z ~ 6 FSRQs would be de¬ 
tectable even at the MeV band. 

5. Discussions and Conclusions 

In this paper, we studied the expected number counts and 
redshift distributions of AGNs (Seyferts and FSRQs) for future 
MeV missions based on recent AGN FFs and spectral models 
as Seyferts (108) and FSRQs (A09) are discussed as the plau¬ 
sible origins of the cosmic MeV gamma-ray background. For 
Seyferts, we assume two primary spectral models. We consid¬ 
ered a thermal spectral model with a cutoff at 0.3 MeV (U03) 
and a thermal plus non-thermal spectral model (108). We adopt 
the (U03) FDDE XFF for both Seyfert models. For the ther¬ 
mal plus non-thermal spectral model (108), a non-thermal com¬ 
ponent appears after the thermal cutoff at 0.3 MeV with the 
photon index of 2.8 to explain the MeV background up to a 
few tens of MeV. For FSRQs, we adopt two models. One is 
based on the Swift -BAT detected FSRQs (A09) and the other 
is based on the Fermi-LAT detected FSRQs (A12). The A09 
FSRQ model can explain the entire MeV background solely by 
FSRQs, while the A12 FSRQ model makes up up to ~ 30% of 
the MeV background by FSRQs. 

Since a thermal cutoff exists in spectra in the U03 Seyfert 
model, we can not expect any detections at > 1 MeV even with 
the sensitivity limit of 10“ 12 erg cm -2 s _1 . In contrast, if the 
origin of the MeV background is non-thermal tail from Seyferts 
(108), we can expect several hundred Seyferts at the MeV band 
with the sensitivity limit of 10” 12 erg cm -2 s -1 . Since each 


Seyfert is faint, we can detect only nearby (z 1) Seyferts even 
with the sensitivity limit of 10 -12 erg cm -2 s -1 . 

If FSRQs make up the whole MeV background (A09), the 
sensitivity of ~ 4 x 10 -12 erg cm -2 s -1 is needed to de¬ 
tect several hundreds of FSRQs at the MeV gamma-ray band. 
However, based on the latest FSRQ GFF (A12), the sensitiv¬ 
ity limit of 10~ n erg cm -2 s -1 , which is almost the same as 
the expected sensitivity of the next generation MeV telescopes, 
would be enough to detect several hundreds of FSRQs. The 
difference between the two FSRQ models comes from the dif¬ 
ferent evolutionary history. 

Future MeV observational windows will range about three 
orders of magnitude in energy. With that wide energy range, 
the dominant population of the cosmic MeV gamma-ray back¬ 
ground radiation would change with the energy. Furthermore, 
as there are uncertainties of luminosity functions and spectral 
models, the sensitivity of several times of 10~ 12 erg cm" 2 s -1 
would be desirable to detect several hundred AGNs. 

It is important to know the origin of the MeV background 
because the expected source counts depends on it. To observa- 
tionally unveil the origin of the MeV background, we need to 
resolve the sky into point sources as in soft X-ray. Here, the ex¬ 
pected sensitivity of the next generation MeV instruments such 
as ASTROGAM is ~ 10~ n erg cm -2 s -1 . With this sensi¬ 
tivity, the expected resolved fraction will be a few percent of 
the total background in the 108 Seyfert scenario, while it will 
be ~5% in the A09 FSRQ scenario. Thus, the large fraction 
of the MeV sky will not be resolved even with the next gener¬ 
ation instruments. However, as about one hundred to several 
hundred FSRQs are expected, we will be able to obtain typical 
MeV spectra and cosmological evolution of FSRQs and the¬ 
oretically estimate their contribution to the MeV background 
more robustly. On the other hand, it would be hard for Seyferts 
since we can expect only about ten sources. 

The angular power spectrum of the MeV sky will provide a 
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unique opportunity to unveil the origin of the MeV sky (Inoue 
et al., 2013). The SGD onboard ASTRO-H will be able to probe 
it up to 0.6 MeV. Although the angular power spectrum with the 
sensitivity of the SGD is dominated by the Poisson term Cf, a 
significant difference of the C? of Seyferts (Inoue et al., 2008) 
and FSRQs (Ajello et al., 2009) is expected. This is because 
FSRQs are brighter but fewer than Seyferts. Anisotropies of 
the MeV sky would allow us to understand the origin of the 
MeV background and be useful to design future missions. 
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ful reading and constructive comments. The authors also thank 
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ments. Y.I. acknowledges support by the JAXA international 
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